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Abstract. We present a detailed analysis of one-dimensional Neel walls in thin 
uniaxial ferromagnetic films in the presence of an in-plane applied external field 
in the direction normal to the easy axis. Within the reduced one-dimensional 
thin film model, we formulate a non-local variational problem whose minimizers 
are given by one-dimensional Neel wall profiles. We prove existence, uniqueness 
(up to translations and refiections), regularity, strict monotonicity and the precise 
asymptotics of the decay of the minimizers in the considered variational problem. 
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1. Introduction 

Ferromagnetic materials are at the heart of modern information storage technology, 
whose need to keep up with the ever-growing amount of digital data, currently in 
excess of 10^^ bytes worldwide, is readily apparent [T]. This is why these materials 
have attracted a huge degree of attention since the early days of the digital age. 
The basic principle of magnetic storage relies on the tendency of the electron spins 
in ferromagnetic materials to align along certain preferred directions, giving rise to 
magnetic domains [2]. Registering and manipulating the magnetization orientation 
in a given domain is then used to read and write discrete data encoded by the 
magnetization orientation in each domain. 

One common magnetic storage solution relies on the use of thin uniaxial 
ferromagnetic films in which the magnetization vector prefers to align along either 
direction of the easy magnetocrystalline axis in the film plane [TH5]. When the 
film thickness becomes sufficiently small (under a few tens of nanometers), the 
magnetization vector is constrained to lie almost entirely in the film plane. In this 
situation magnetic domains in epitaxial (monocrystalline) films usually consist of 
relatively large regions, in which the magnetization vector is nearly constant and 
oriented in the direction of one of the two possible directions along the easy axis. These 
regions are separated by narrow transition regions, called domain walls, in which the 
magnetization vector rapidly rotates between the two orientations [HISHS] • One of the 
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most common wall types in such materials is the Neel wall, which separates two regions 
of opposite magnetization by an in-plane rotation and is oriented along the easy axis 
to ensure zero net magnetic charge. In real materials these walls are often pinned to 
the material imperfections, and their motion determines magnetization reversal under 
the action of applied magnetic fields [2] . 

Studies of Neel walls have a long and somewhat controversial history (see the 
discussions in [2j[T0]), but at present the structure of the Neel wall in very thin 
films appears to be rather well understood on the basis of micromagnetic arguments 
PIMTTHTT] . The basic features of the predicted one-dimensional Neel wall profiles had 
been verified experimentally [18]. Rigorous mathematical studies of the Neel walls 
are more recent and go back to the work of Garcia- Cervera [T4l[T6], who undertook 
some analysis of the associated one-dimensional variational problems and performed 
extensive numerical studies of the energy functional obtained by Aharoni from the full 
micromagnetic energy after restricting the admissible configurations to profiles which 
depend only on one spatial variable |19| . Melcher further studied the minimizers 
of the same functional in the class of magnetization configurations constrained to 
the film plane and established symmetry and monotonicity of the energy minimizing 
profiles connecting the two opposite directions of the easy axis [20]. Using a further 
one-dimensional thin film reduction of the micromagnetic energy introduced in |17| . 
Capella, Melcher and Otto outlined the proof of uniqueness of the Neel wall profile and 
its linearized stability with respect to one-dimensional perturbations [21j . Stability 
of geometrically constrained one-dimensional Neel walls with respect to large two- 
dimensional perturbations in soft materials was demonstrated asymptotically in |22j . 
More recently, F-convergence studies of the one-dimensional wall energy in the limit of 
very soft films and in the presence of an applied in-plane field normal to the easy axis 
were undertaken in [23l[24], and a rigorous derivation of the effective magnetization 
dynamics driven by the reduced thin film energy introduced in |21j from the full 
three-dimensional Landau-Lifshitz-Gilbert equation was presented in [25] . 

In this paper, we perform a detailed variational study of the Neel walls, 
understood as one-dimensional minimizers of the reduced thin film micromagnetic 
energy, in uniaxial materials in the presence of an applied in-plane magnetic field in 
the direction perpendicular to the easy axis, extending previous results for Neel walls 
in the absence of the applied field. We prove existence, uniqueness (up to translations 
and reflections), regularity, strict monotonicity and the precise decay of the energy 
minimizing wall profiles. Our variational setting is slightly different from that adopted 
in the earlier works and relies on the angle variable rather than the two-dimensional 
unit vector representation of the magnetization. For this reason our proofs differ in 
a few technical aspects from those of [20]. In fact, one of the purposes of our work 
was to clarify some of the arguments in the analyses of [16l[20l[2T] . In particular, we 
spell out the details of the proof of uniqueness of minimizers within our setting and fill 
in the missing argument for proving strict monotonicity of the angle variable as the 
function of coordinate, which is needed to establish stability of the Neel wall profile 
in [21] . We also establish the precise asymptotic behavior of the Neel wall profiles at 
infinity, which is new even in the case of zero applied field. Let us note that while in 
this paper we are not concerned with the logarithmic tail of the Neel walls in very soft 
materials, which was one of the main focuses of [TU[TUIin], our decay estimates could 
be made quantitative in this regime to yield the intermediate asymptotics of the Neel 
wall profile away from the core. 

The rest of our paper is organized as follows. In Sec. [2] we discuss the basic 
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micromagnetic energy and derive the reduced one-dimensional energy that describes 
the Neel walls in the applied in-plane field oriented normally to the easy axis. Then 
in Sec. [3] we present the variational setting for our analysis and state our main result. 
Sec. 2] contains a few auxiliary results and Sec. [5] contains the proof of the main 
theorem. We also discuss some open problems at the very end of Sec. [5l 

2. Model 

In this paper we are interested in the analysis of the energy minimizing magnetization 
configurations in thin uniaxial ferromagnetic films of large extent with the easy axis 
in the film plane. We also wish to include the effect of an applied in-plane field 
in the direction normal to the easy axis. The starting point in the studies of such 
systems is the energy functional, introduced by Landau and Lifshitz, which leads to a 
non-convex, nonlocal variational problem. The functional, written in the CGS units, 
is [3[ini[ll[27]: 

Here f7 C M"^ is the domain occupied by the ferromagnetic material, M : M'^ is 
the magnetization vector that satisfies |M| = Mg in Q, and M = in IR'^\ri, the positive 
constants M^, A and K are the material parameters referred to as the saturation 
magnetization, exchange constant and the anisotropy constant, respectively, Hcxt is 
an applied external field, and $ : R'^ — R is a non- negative potential that has several 
minima at which $ vanishes. Note that V • M in the double integral is understood in 
the distributional sense. 

The micromagnetic energy in ([T]) is composed of five terms: the exchange energy 
term which penalizes the spacial variations of the magnetization M, the anisotropy 
term reflecting the magnetocrystalline properties of the material, the Zeeman energy 
favoring the alignment of M with the applied external field, the stray-field energy, 
which is nonlocal and favors vanishing distributional divergence, i.e., V • M = in 
both Q, and d^l (the so-called pole-avoidance principle), and an inessential constant 
term added for convenience. In the case of a uniaxial material of interest to us, 
there exists a distinguished axis identified through a unit vector e, and $ is given by 
$(M) = - (M • e)2, so that the minima of $ are {ieMg} pi llOl[?7] . 

In the case of extended monocrystallinc thin films with the in-plane easy axis we 
have J7 = X (0,c?), and without loss of generality we may assume that e = 62, 
where is the unit vector in the i-th coordinate direction. For thin films (moderately 
soft, ultra-thin) of practical interest to magnetic device applications such as MRAMs 
(magnetoresistive random access memories) !4-6 , a significant reduction of the energy 
in ([l} is possible, giving rise to the reduced thin film energy [9lll7j . To better explain 
the relevant parameter regime, let us introduce the following quantities 



called the exchange length, the Bloch wall thickness, and the material quality factor, 
respectively. When the film is ultra-thin and soft, we have d < £ L, but at the same 
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time one also has a balance Ldjf? ^ 1 for many materials [5H]. In this situation the 
dimensionless parameter 

_ 47rAf2(i _ Ld _ d 

which is referred to as the thin film parameter [17j , becomes a single measure of the 
strength of the magnetostatic interaction relative to both anisotropy and exchange. 

The reduced thin film energy is formally obtained from the full micromagnetic 
energy in ^ by assuming that M does not vary in the direction of 63 (the direction 
normal to the film) , setting the component of M along 63 to zero and passing to the 
limit Q — >■ and d — ?> jointly, subject to = 0{1) fixed, after rescaling lengths 
with L [17] (see also [9]). Assuming further that Hcxt = eihK/Ms, after a suitable 
rescaling we arrive at the following reduced energy functional: 

E{m) = ]-( \\/m\^d^r+l [ {m-ei-hfd^r 
2 Jr2 2 Jjj2 

f /" V ■ m(i")V • m(r') 



in 



d^rd^r', (4) 



where now m : — > §^ is the unit vector in the direction of the magnetization in the 
film plane. Note that the assumptions on M used in this derivation are justified by 
the strong penalization of the variations of M across the normal direction to the film 
by the exchange energy and by the strong penalization of the normal component of 
M by the shape anisotropy [2]- Also, up to a constant factor the last term in (j4]) is 
simply the square of the homogeneous 7J~^/^-norm of V • m in [9]. 

The reduced energy in ^ is the starting point of the analysis of the rest of our 
paper. Without loss of generality we may assume that h> 0. Also note that for /i > 1 
the energy in (|4]) admits only a unique global minimum m = ei and no Neel walls are, 
therefore, possible in this situation. For h € [0, 1), on the other hand, there are two 
global minimizers m-|- = (/i, ±-\/l — /i^) corresponding to the two monodomain states. 
In the following, we will always assume that h is in this non-trivial range, in which 
Neel walls connecting the two states appear. Let us point out that at the same time 
we do not allow the external field to have a component in the direction of the easy 
axis, since in this case only one monodomain state exists as the global minimizer of 
the energy. Under an applied field in the direction of the easy axis Neel walls begin to 
move, invading the domain with higher energy density by the domain with the lower 
energy density [HHI]. Similarly, the considered wall orientation along the easy axis 
is the only one that makes the stray field energy of a one-dimensional profiles finite. 
When the wall makes a non-zero angle with the easy axis (cf. with [29]), it carries a 
net magnetic charge, which makes the associated magnetostatic potential for the wall 
in the whole of infinite. 



3. Variational formulation and statement of the main result 

We now turn to the study of one-dimensional Neel wall profiles. For that we assume 
that m varies only along ei and compute the energy of such a configuration per unit 
length of the wall. It is convenient to introduce the new variable ■& = ^{x) which 
gives the angle that the vector m makes with 62 in the counter-clockwise direction as 
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a function of the coordinate along ei. Thus, setting 

m(x) = (-sin?9(x),cosi?(a;)) e §^ (5) 

for every a; e M, we can rewrite the one-dimensional version of the functional in Q in 
terms of the angle variable d to obtain the one-dimensional Neel wall energy (cf. [ITj): 

E{^;R):^^J ||^^|2 + (sin?9 - /i)^ + ^ sini?(^ - ^ sini^jda; 



, 2 



(6) 

1 /2 

where, as usual, (— d^/dx^) denotes the square root of the one-dimensional negative 
Laplacian (a linear operator whose Fourier symbol is |fc|), and, furthermore, we used 
the identity [29 1 130 ] 

. ^ 1 i u{x) - u{y) 

— u{x) = -f— T^dy, (7) 

dx ) ttJr [x — yY 

for every x and, say, every u S C^(R), where _f stands for the principal value of the 
integral. 

We wish to study the minimizers of the energy in (|6]) among the profiles that 
connect the two distinct minima of the energy at x = ±cx). To this end, we need 
to introduce a suitable admissible class of functions which yields minimizers with the 
desired properties. We propose to minimize E{;d\ K) over the admissible class 

A:^{d^ulM)--'^-'nh^ii^m\, (8) 

where r\h G C°°(M; [0, 1]) is a fixed non-increasing function such that, setting 

61/1 ~arcsin/iG [0, f), (9) 

we have rjh = tt — 0h in (— oo,— 1) and rjh — 6^ in (l,+oo). We point out that the 
definition of does not depend on the choice oirjh'- If G C°°(R; [0, 1]) is a different 
non-increasing function, i.e., f]h 7^ satisfying T)h — — Oh 'va. (—00, —1) and T)h = Oh 
in (1, +00), then 

A^{de i/L W --^-flh^ H\R)}. (10) 

Indeed, any d d A satisfies also — fjh E H^{R) for ry/j — rjh E H^{R); vice versa, 
for the same reason any 1} E H^-^^CR) with d — fjh E H^{K) belongs to A. Note that 
our choice of the admissible class A fixes the rotation sense of the Neel wall, and the 
wall of the opposite rotation sense may be obtained from the minimizer over ^ by a 
reflection about a; = 0. 

It is easy to see that the Euler-Lagrange equation associated with the functional 
in (|6]) is given by 



V „ / d 



•2 \ 1/2 
2 V dx^ 



cos 1? sin I? — /i cos 1? + — cos 1? ( 3— ) sin?? = 0, (11) 
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with the boundary conditions at infinity 

lim d{x)^eh, lim d{x)='K-eh. (12) 

x— f+oo a:— > — oo 

The main result of this paper is the following. 

Theorem 1 (existence, uniqueness, regularity, strict monotonicity and 
decay of Neel walls). For every u > Q and every /i G [0, 1) there exists a minimizer 
of E{'d;M.) in ^ over A in fS]), which is unique (up to translations), strictly decreasing 
with range equal to {9h,'^ — 0h) and is a smooth solution of that satisfies the limit 
conditions given in (jl2[) . Moreover, if : R — ^ (9h, tt — 9h) is the minimizer of E 
in the class A satisfying '[^(''^(O) = ^, then •d'^'^^{x) = it — 'd^'^\—x), and there exists a 
constant c > such that lim^^-^+oo x"^ {i)''^^ {x) — Oh) = c. 

4. Some auxiliary lemmas 

We start with a few preliminary considerations and lemmas. Let -d ^ A. By Morrey's 
Theorem (see [3T1 Theorem 11.34]), i9 - m ^ C^^'^{M.) and - 7/^ ^ as x -> ±oo; 
that is, 7? S C(R) H L°°(R) and satisfies ((HI). Furthermore, assuming in addition that 
i9(R) C [6lft,7r - 0h] and E{^,R) < +oo, and defining p:R~^ [Oh, f ] by 

ifi?(:.)e[0ft,f], 

for every a; G M, we have sin p — sini9, and since the map i9 i— > p is Lipschitz, we also 
have = \Px\ almost everywhere on M. Thus 

+0O > Ei'd, M) = E{p, ^) = IJ^ Ip-I" dx + ^ Jj^'^'^'P' dx + j\\ sinpll^^,,,^^^ 



-2 jj^''^ '^^^T in -||sinp||^i/,^^^ 

^ x"^"^''"^^! + ^^^^ 

where c/i := cos^ (f + > for all < ^. In the above inequality, we used the 
fact that, since p(M) C [0h, f]. Oh G [0, f ), and sinz > z/V2 for all z e [0, f ], so that 

sinp - h = sinp - sin6',i = 2cos (^^~~^^^ (^~2~^) ~ v^'''^ ~ ^^^^ 

Lemma 2 (restriction of rotations). Let -d A such that E(i3) < +oo. Then there 
exists '& G A such that ■d{R) C [6'/j,7r — Oh] and E{'d) < E{'d), with strict inequality 
unless 7?(M) C [Oh,TT - Oh]. 

Proof Step 1. We show first that there exists 'd^ £ A such that ■d'^{R) C [0,7r], 
sin^?^ |sin??|, and E{d'') < E{d). Let : R [0,tt] be defined by 

, , r ^{x)-2kn if 7?(x) e [2fc^, (2fc + 1)^) 

\ 2kTT-i){x) if 7?(a;) e [(2fc- l)7r,2fc7r), ^ 

for every x e M. The definition is well-posed since { [2fc7r, (2fc + l)7r), [{2k — l)7r, 2fc7r) : 
/c € Z} is a partition of IR.. Notice that -i?"^ is obtained by means of a translation of 
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—2kTT in each interval of the form [2fc7r, (2fc + l)7r), and by means of a reflection 
with respect to the origin and a translation of 2fc7r in each interval of the form 
[{2k — l)7r, 2fc7r). By construction, ■d'^ G [Oj'"'] and sin?? < |sin??| = sin??'^ so 
that II sin??||i2(R') = || sin'!9'^||i2(][j') and —h J^shi'd'^dx < —h J^sinddx, implying that 
J^{sm-d'^ — h)'^dx < J^{sm'd — h)'^dx. Furthermore, since the map i— i?"^ is Lipschitz, 

we also have J^-d^dx = (z?^)^da;. 

The conclusion now comes from the fact that 

d' V^' , 1 f f {u{x)-u{y))\ 
' udx — — f f — ax dy 



dx^ J 27r 7r 7r (a; - y)2 

JmJr [x — y) Jr \ ax / 

Step 2. Without loss of generality, in view of Step 1, we can assume that d £ A is 
such that t?(R) C [0,7r] and £;(i?) < +oo. Set / {a; e R : ^{x) e [0, 6*^) U (7r-6l,,, vr]} 
and let d: M. ^ [9h,n — 9h] be defined by 

r Oh \i'&{x)<eh 
^{x):=\ d{x) ii§{x)e[dh,Ti-eh]=^\I, (18) 

{ ^-Oh if ^9(a;) > ^ - dh 

for every x e M. 

Notice that, since i? = i? in R \ /, we have 

II sin^ - h\\l.^^~^ - II sin,? - h\\l.^^) = II sin^ - h\\l,^,^ > 0. (19) 
Moreover, since i? is constant on /, we have 

If f ( sin 19(2;) — sin ??(y))^ 



w^^^h.^m - 2^4,4, — — '^'y 

If f (h- sm-d(y))^ , ^ 



and so 



„2 „ ~,,2 1 f /■ (sinz9(a;) - sini?(y))^ 
I|sinz9|l^,/,(j,)-||smi?||^,/,(jj^ = — jj^ ^— dxdy 

+ i/ /(sin^(.)-sin%))- (.-sin^(,))-^^^^ 

I" Js.\i J I - yY 

1 f /■ (sin??(x) - /i)(/i + sin??(x) - 2sinz9(y)) 
> - / / ^ '-^ "2 dx dy > 0, 

where in the second to last inequality we have used the identity A^ — B^ = {A — B){A+ 
B) and in the last one the fact that for every x € I and every y e R \ / the following 
inequalities hold: sin 79(0:) < h < sm'0{y). Then we conclude that < E{d). □ 

The following rearrangement property is a consequence of Lemma 7.17 in |32) . 
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Lemma 3. Let u G L^(R) he a nonnegative function and let u* he its symmetric 
decreasing rearrangement, i.e., 

u*{x) / X{u>tr ix)dt, (22) 
Jo 

where for a Borel set A C M the rearranged set A* is the interval with measure C^{A) 
centered at the origin. Then 

d^ _ ^.^2x1/2 



ll"*lll,i/.(K) u*dx< J^ui^-^j udx=\\u\\%,,^^y (23) 

with equality only if u is a translation of a symmetric decreasing function. 

To prove our main result, we first need the following preliminary lemma. The idea 
goes back to [20], except that here our variable is the angle function rather than the 
first component of the magnetization vector. Let us point out that the rearrangement 
argument of Lemma |4] only yields the non-increasing property of 'd for the minimizer. 
To prove strict decrease, one needs an additional argument presented in Step 3 of the 
proof of the main theorem below. 

Lemma 4 (rearrangement). Let E A be such that z9(M) C [Oh,TT — Oh] and 
E^t}) < +O0. Then there exists a function 'd°{x): M — > [^?/i,7r — 9h] satisfying (|12p 
and the following properties: 

r(0) = f , i)°{x) ^TT-d°i~x), i9°<0 on M and Ei-d") < E{i}), (24) 

where the equality in the latter expression holds only if sini? is a translation of a 
symmetric decreasing function. 

Proof. Let p: M. —i' [6h, ^] he defined as in ((T5| for every a; € M. Then, from the 
discussion at the beginning of Sec. [5] we have E{d,R) = E{p,'M.). Now, define 
p° : R [6*,,, f ] by setting 

p°{x):=9h + {p(x)-9hy, (25) 

where given a function /, /* stands for the symmetric rearrangement of /. This 
implies that p° is even, {p*)x < on K+, and p°{x) — >■ 9h as |a;| — > +oo. Moreover, 
the level sets of p° are simply the rearrangement of the level sets of p, i.e., 

{x: p''{x)>t} = {x: p{x)>t}*. (26) 

A consequence of this is the equimeasurability of the functions p° and p, i.e., 

C^{{x : p°{x) > t}) = : p{x) > t}) (27) 

for every t > 0. This, together with the Layer Cake Representation Theorem 1.13 
in [32], yields 

/ (j){p''{x))dx= [ (j){p{x))dx (28) 
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for every monotone, absolutely continuous function <j) : [0,00) — )■ [0,00) satisfying 
(j){0) = 0. Choosing (f>{z) — [z — h)\, we see that 



{^ sva p° {x) — h^"^ dx — i ( sinp — /i)^(ia; 



(29) 



We recall now the following property of rearrangements of any given Borel 
measurable function / : R R vanishing at infinity: 



($01/1)* =$0/* 



(30) 



for every $ : R+ — > R+ nondecreasing. Applying the above property to / := p — Oh ■ 
M — >■ [0, ^ — Oil] and sin{9h + /), which is increasing for every / G [0, f — 9h], 

we get 

(sinp)* = sm{9h + [p - 9^)*) = smp°. (31) 
In view of the above identity and using also Lemma [3l we have 

f f d'^ V^'^ f f <P 

j sinp°( --^1 sinp°da:: = j (sinp)*( --^j (sin p)* da; 



1/2 



siB.'d dx. 



(32) 



where the equality holds only if sinz? is a translation of a symmetric decreasing 
function. Finally, by Lemma 7.17 in [32], 



ol?dx = 



{p(x) - 9nyX dx < I / 1^9, 



'dx. 



Define : R ^ [6*^, vr - 9^] by setting 
d°{x) :-- 



p°ix) 



if x > 
if a; < 0. 



Since smd° = smp°, by ((29]), ([32]), and ((33]) we conclude. 



(33) 



(34) 



□ 



We now investigate the decay of monotone solutions of ([TT|) satisfying (fT2|). 
This information, combined with the properties of the fundamental solution of the 
linearization of (fTTj) around 9h (see, e.g., [Ml Sec. 5.1]), will be used to establish the 
precise asymptotics behavior of the minimizers of E over ^ as a; — >■ ±00. 

Lemma 5. Let i9 e C°°(R) be a non-increasing solution of (|lip satisfying (jl2p . wit/i 
'(9(— a;) = vr — ?9(x) /or x G R and all derivatives vanishing as x ^ ±00. Let 
u = sin?? — h and assume that there exist c > and a G (0, 2] sitc/i i/iat 



M(a;) < 



Va; G 



(35) 



T/ien there exists C = C{c, a, h,iy) > such that 



\Uxix)\, \Uxxix)\, 



dx"^ 



1/2 



u{x) 



< 



_C 

1 + kl" 



Va; G 



(36) 
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Proof. Throughout the proof we assume that x > is sufficiently large depending only 
on c, a and h. All the constants in the estimates are also assumed to depend only on 
c, a, h and v. 

Equation (|lip written in terms of u reads 



{h + u)ul 



l-h2-2uh-u^ (1 - /i2 _ 2u/i - 1*2)2 



(ix2 



1/2 



(37) 



In particular, since u{x) goes to zero as a; — >■ oo together with all its derivatives, we 
also have that lim M{x) — 0, where 



M{x) 



max 

[x,oc) 



1/2 

dx"^ 



(38) 



Now, multiply ()lip by "dx and integrate over (x,oo). Together with the monotonicity 
of u(x), this yields 



-dlix) < u'^ix) + iyu{x)M{x). 
Using the fact that Ux — cos from ((37)) and (p9| we obtain 

\uxx{x)\ < C{u{x) + M{x)). 
On the other hand, for every S > sufficiently small we have by ^ 



(39) 
(40) 



d2 



1/2 



'd^ ' ^ - 



x — 6 



(x-7y)2 



"2^ u(a;) - u{y) ^ 1_ 
{x- yY TT 

u{x) ~ u{y) 
x~s {x - y)2 



2^ u{x) - u{y) 

—'^y 



dy+- 



1 

2^ 

+ 00 



x+5 



{x - yf 

u{x) - u{y) 
{x - yY 



dy. (41) 



Using the symmetric decreasing property of u{x)^ we can then estimate the left-hand 
side of (HH) as 



da; 2 



u{x) 



<C\6 ^u{^x) + 5 max + 

[x — 5.x+5\ 



\x {x-yY 



dy 



(42) 



where to estimate the fourth term in the reight-hand side of (1411) we used Taylor 
formula u(y) — u{x) + Ux{x){y — x) + ^Uxx{x{y)){y — xY , for some i(?;) lying between 
X and y, and noted that the linear term in y — x does not contribute to the principal 
value of the integral. Combining the estimate in (j42|) with (j40| and the assumption 
on the decay of u in ([SS]) then yields 



da; 2 



1/2 



u{x) 



< C 



— + 6Mix -6) + — + 



l+a 



Now, choosing S > sufficiently small, the last estimate implies 

d2 X 1/2 



da;2 



M(a;) 



< Cx'" + ^M{x - 5), 



(43) 



(44) 
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for every a < 2. Therefore, by monotonicity of M{x), we also have 

M{x) < Cx-" + \M{x - (5), (45) 

for all X > sufficiently large. 

Let us show that implies the same kind of bound as in ([55)1 for M{x). 
We use an induction argument. Let x„ := xq + n5 and Af„ := M{xn) for n e N 
and some xq > to be fixed shortly. Clearly, by ps)) we have Mi < cx^°' for 
c = 4C + iAf(a;o)(a;o + 6)°". We claim that if also M„_i < ca::^"i, then (gS]) implies 
M„ < ca;~", provided that xq is chosen to be sufficiently large. Indeed, by dJS]) and 
the assumption Mn-i < cx~"i we have 

Mn<(c+ ^(^^^^) 2;-" < (C + |c)x;;" < (46) 

provided that xq > 66, and the claim follows. Once we established the bound for 
Mn, the estimate M{x) < C/(l + x") for all a; e M follows by monotonicity and 
boundedness of M{x). 

To conclude the proof of the lemma, we combine the estimate for AI[x) just 
obtained with ^ and dM]). □ 



5. Proof of the main result 

We are now in a position to prove the main result of this paper. Theorem [TJ 
Step 1: existence. For every n let C .4 be a minimizing sequence, i.e., 

lim E{^n) = M{E{d) :^£A}< +oo. (47) 

By translation invariance we may assume that (0) = ^ and by Lemma S] that each 
■dn satisfies 

M^) ^[eh,7T-9i,], dn{x) ^TT-dn{-x), lim l9„(x)=0ft, (l?„), < On M. 

a: —^+00 

(48) 

For each n, let p„ : M — > [0, ^] be defined as in (fT3|) (after replacing i) with i9„). By 
(fn|) we have 

- 6'/i||hi(m) + II sinp„|||i/2(R) <C < +00. (49) 

In view of Banach-Alaoglu-Bourbaki's Theorem, there exist a subsequence of {p„}, 
not relabeled, a function p € i/^(R) and a function u S i?^/^(R) such that 

Pn — f^h ^ V weakly in H^{W) and sinp„ ^ u weakly in _ff^/^(R). (50) 

Let us fix fc e N. By the compact embedding H^{—k, k) CC L'^{~k, k), we may 
find a subsequence, not relabeled, such that pn — Oh ^ v strongly in L'^{—k,k) and 
£-'^-a.e. in {—k,k). Thus, p„ — > p £^-a.c. in {—k,k), where p := v + Oh- Moreover, 
since 

sup||smp„|!^,,,^_^^^^=sup(^- y_J_^ ^— p dxdyj < +00, 

(51) 
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it follows by the Fractional Compact Embedding (see for instance Equation (5) in 
Section 8.1 of [32]) that {sinp„ — h} is precompact in L'^{—k,k), that is, up to a 
subsequence, not relabeled, sinp„ u strongly in L'^{—k,k) and £^-a.e. in {—k,k). 
Then by the uniqueness of the limits, u — smp C^-a.e. in {—k,k) for every k. 

Finally, by the lower semicontinuity of the and i/^/^ norms with respect to 
their weak convergences and Fatou Lemma applied to the sequence {(sinp„ — h)^}, 
we have 

1 /•', ,2,^1 f% ■ ,,2 , ^ ^ f' t (sinp(x)-sinp(j/))^ ^ 

<liminf- / \{pn)x\'^ dx + - / (sin p„ - /i)^ da; (52) 

n^+oo 2 J_k 2 

+f r r ^""^"^;^'^;7"^'^^' ^^^^<""^infi^(p„,R). 

Applying Lebesgue's Monotone Convergence Theorem to the sequences {x{~k,k) [(^2;)^+ 
(sinp- /i)2]} and |x(-fc,fc)" ^^'"''^('il^')"'^^^^^ }' '^^ ^'^'^'^ obtain 

£;(p,M) < liminf i;(p„,R) = liminf £;(i^„,R). (53) 

n— ^+00 n— J-+00 

Given such a p, let us define : E — > M by setting 

^ ^ I ^ - ^(2;) if X < 0, ^ ^ 

for every x G M. We claim that satisfies the following properties: 
^)(o)(0) = f, i?(o)(a;) =7r-?9(o)(-a;), lim (x) = 6*^, < on M. (55) 

a;— >+oo 

Since p — 9h G 77^ (R), by Morrey's Theorem we have lim 'd^^'^{x) = lim p{x) — 9h 

and e C(R). Finally, since weak convergence in H^{M.) implies pointwise and 
uniform convergence on compacts (see for instance Theorem 8.6 in |32j). by the first 
of (f50|) we have pn ^ p locally uniformly on R. Therefore, the following properties 

Pn(0) = f , pn{x) ^ Pn{-X), Pn{xi) > Pn{x2) if ^ < Xi < X2, (56) 

are preserved when taking the limit as n — ^ +00 which in turn implies by construction 
that the properties in (l55t are satisfied. 

Since by construction i?(i?(°),R) = E{p,M), by and ([51]), we conclude that 
f^'^) is a minimizer for E in the class A. 

Step 2: Euler-Lagrange equation. Since the first variation of E at any global 
minimizer of E is zero, we have that the minimizer i?'"-* of E^ as well as any other 
minimizer of E, satisfies 



dt 



+ (f COS 

^(0) sin,9(") - /i(pcost?(°) 

+ —ip COS 

^(0)(^_^^ sin^(0)^da;, (57) 
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for every Lp £ iJ^(R). In other words, (f) = is a weak solution of the ordinary 
differential equation 

— (fir^x + b{x) COS = 0, (58) 

where b(x) := sini?(o) -h + ^v{-d^ /dx^)^/^ sm'd^°'> . Since S L^iR) and the sine 
function is Lipschitz, we have (sini?'^"))^: = cosi5*^°^ G L^(M), and by the estimate 

<C\\u\\Hum, (59) 



/ j2 \ 1/2 



2 V da; 



L2( 



applied to u := sini?^'^^ — h, we conclude that 6(a;) € L^(M). By (|57|) . the weak second 
derivative of is given by -d^xx — b{x) cos 1)^^'' and is, hence, in L^(]R). Therefore, 
■&^x ^ e iJi(R) and by Morrey's Theorem, together with the fact that i?^") G C(R), 
we have that "dl^"^ is continuous and bounded on R as well. Using the interpolation 
inequality (here we follow the arguments of [21 ) 

W^^^^hnm < ll'Ji^^lli^iwIl^i^^lliiw < +°«' (^o) 

we have (simJ^"))^^ ^ iji^j cos i?(o) - (iJI^^ sini?(°) E L^{M.). Thus, by ^ and 
the fact that {d/ dx){—d'^ / dx'^Y^'^u = (— (i^/dx^)^/^Wa;, readily verified via the Fourier 
representation of the fractional Laplacian [30], we conclude that b{x) € H^{M) and 
so by Morrey's Theorem b{x) is continuous and bounded. This in turn implies 
that d^xl = -6(x)cosi?(°) e L°°(R) n C(M), that is, z?'"' e C2(R) is a classical 
solution of ([58]) . Finally, bootstrapping these regularity arguments, we conclude 
that e C°°(R), and, moreover, all the derivatives are bounded uniformly in M. 
Together with boundedness of their L'^ norms, the latter implies that all the derivatives 
of i?'^"-' vanish at infinity. 

Step 3: strict monotonicity. First we show that d^x\^) < employing the same 
argument as in [21]. Since i^^^-* G C^{R) is a classical solution of ([SS]) . if we assume 
that ■(?(°)(0) = f and ^^x\o) = 0, the uniqueness theorem for ordinary differential 
equations implies that 1}^°^ must be identically equal to |, a contradiction. 

Now we prove /x ^ < OonR+. We argue by contradiction and assume that there 
exist a; > such that -dl^^x) = 0. But then also -d^xx must be zero at x, because 
otherwise 1}^^^ will be either strictly convex or strictly concave at x, contradicting the 
fact that 'd^x'^ < on R. Moreover, by differentiating the Euler-Lagrange equation 
(|lip and taking into account that ^x{x) — ^xx{x) = 0, we get 

eU^;) = ^cos/0)(i)g(x), g := - ^^^^ (^o) ^(0)), (gl) 

Now we observe that cosz9(°)(x) > because ??(°)(R+) C [9h, f ], and, recalling (O, 
taking into account that "dxix) = ^xxix) — 0, and noticing that the integral converges 
near x in the usual sense, we have 



9{x) = / Tz 72 dy 



1 /■ ^i°)(2/)cosi?(0)(y) 1 [ ^^\y)cosm;yl 

■dy - - I , ,2 dy > 0, (62) 



TT Jr+ {x - yY TT {X + y) 
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where we have used the fact that 'd^x\—x) = 'd'"x\x) < 0, and the inequahty is strict 
on a set of positive measure. Therefore, "dl^Jxix) > 0, which implies that i?'") is 
increasing at x, contradicting the fact that ■iji"'' < on M. Hence, 1?'-°' has to be 
strictly decreasing on R+. 

Since d'-^^x) = n - for every x < 0, it follows that z?!"^ < also on R". 

Step 4: uniqueness (up to translations). It is clear that in view of the 
translational invariance the function = ??(°)(a; — xq) with any xq € R still 

belongs to A and satisfies ?9*^^''^(xo) — ^. To conclude, we have to show that every 
minimizers of i? in is of the form -d^^^x — xq) for some Xq G M. Our argument is 
related to the strict convexity of the integrand in ^ written as a function of u = sin d 
and its derivative noted in [21_. 

We employ the strict monotonicity of minimizers, which implies that for every 
minimizer there is a unique point at which -d — ^. Let and t^^^) \)q two 
different minimizers, which, after a suitable translation, satisfy '(^^"'^•'(O) = ^ and 
t?(2)(0) = |. Define ^(x) := arcsin (^SlRl^lllEl+siR^}^^ for x > and ^{x) := 
IT — arcsin ^212-!^ — (a:)+sini9 — (£)^ x < (the function sini9 is symmetric decreasing 
by Step 1). We claim that for all x 7^ we have 



^'^"^ 4-(sin^(i)+sini?(2))2^ ^2^ ^ ' ^^^^ 

Once the claim is proved, we get that G L2(R), and, moreover, that since the 
anisotropy and the stray-field terms in the energy are quadratic in sini5, we get 
E{'d,R) < i[£;(i?(i),R) + E{'d^'^\M.)], which contradicts the minimality of i^^^^ and 
t9(2). 

Let us come to the proof of For simplicity, we set A :— z?!^-* < 0, 

B := "d^x^ < (by the strict monotonicity proved in Step 3), X cos and 
Y := cost9(2). Noticing that < sin??(*) = Vl - cos2 and that -1 < -Vl - h? < 



C0S,9(1) + COSl?(2))2 (^(1))2 ^ 



(2)N 



X,Y < Vl - h'^ < 1, (|M1) can be written as 

^2^2 + 52^2 + 2ABXY _ {AX + BY 



2 



If X 7^ then, due to the strict monotonicity, (X, Y) 7^ (0, 0) and so the 
denominator in the previous expression is strictly positive. After a few algebraic 
manipulations, ()64p is equivalent to 



f{X,Y) := {B^-A'^){X^^Y^)-AABXY + 2{A'^ + B'^)(l-y^{l-X^){l-Y^)^ > 0. 

(65) 

Setting C := fsq^, D aA'b^ ' ''^^'^ noticing that D > and C'^ = I - D'^, we 



obtain 

fiX,Y) 



A^ + B"^ 



= CiX^ - _ 2XY\/l - C2 + 2 - 2v/(l -X2)(l -y2) 



= [XVTTC -YVT~C) + (y/l-X^ - y/l-YA >0 



(66) 
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Step 5: decay. We claim that by the resuhs of the previous steps the unique 
minimizer z?'-'^' oi E va A satisfies the assumptions of Lemma [S] with a = i. Indeed, 

^(0) 

is a smooth decreasing solution of pT|) satisfying ((T2|) with all derivatives 
vanishing at infinity and obeying the required symmetry property. Furthermore, since 
u = sini?'^°^ — /i e L^(]R) is symmetric decreasing, by an elementary property of 
monotone functions (see, e.g., [331 Lemma A. IV]) we have that u{x) < \x\~^/'^\\u\\L2^^y 
Therefore, the conclusions of Lemma[5]apply to We now claim that this, in turn, 
implies the same kind of estimates for p{x) — 6h, where p is defined by (fT3)) with 
§ = Indeed, since = i?i°'cosi?W and u^^ = i^i^i cos - sm'0^°\ 

the estimates for the derivatives follow from (|36l) . and to obtain the estimate for 
{—d^/dx'^Y^'^p, one can use (|42l) with u replaced by p. 
We now rewrite (jll[) in the following form: 

L{p(x)-eu)^ fix), f{x):=f^{x)+f2{x) + f^(x), (67) 

where 

L:=-^ + -.cos^^.(^-— j +cos^^. (68) 

is a linear operator that can be viewed as a map from iJ^(R) to L-^(M), and 

fi{x) := cos9h{cos9h - cos p{x)){p{x) - Oh) 

+ cos p{x){cos 6hip{x) — 9h) — sinp{x) + sinO/i), (69) 

d' X V2 



f2{x) := -vcosOh {--^ ] {cos0h{pix) - 0h) - smp{x) +smeh), (70) 



-12 \ 1/2 



1 / 

/3(a;) 2*^*^'^°^^''^ ^°^'^(^''M ) (sinp(x) - sin6'/i). (71) 

Note that L represents the operator that generates the linearization of (|lip around 
1} = 6h, and all the terms in the definition of / are of "quadratic order" in p — 0^. 
Also note that the fundamental solution G{x) associated with L, i.e., the solution 
of LG{x) = 5[x), where 5{x) is the Dirac delta- function, is a positive, continuous, 
piecewise-smooth function with a jump of the derivative at the origin and a decay 
G{x) ~ \x\~'^ at infinity (see Lemma lA.ip . Moreover, L is invertible, and we have 

p-eh = L-'f{x) = I G{x- y)f{y)dy. (72) 



Observe that since a priori u{x) < c/(l + |x|^/^), using Taylor expansion 
in p — Oh we have |/i(a;)| < C/(l + |a;|^/^)^, and by Lemma [5] we also have 
1/3(2;)! < C/(l + |x|i/2)2. To prove that |/2(a;)| < C/(l + \x\^/^f as weh, we use 
the estimate in (|42l) once again, taking into account that \{pcos6h — sv[\p)xx\ < 
I cos9h — cosp\ \pxx\ + sinp < C/(l + |a;|^/^)^ for x 7^ and arguing as in Lemma 
[5] Thus, we have |/(a;)| < C/{1 + |a;|), and hence using ((72)) . we obtain 
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In view of (|73p . we have now improved the estimate for u{x) to the one in (j35p with 
a = 1. Repeating the argument above, we then conclude that u[x) < C"/(l + 
Let us show that this estimate is, in fact, optimal, and obtain the precise asymptotics 
of the decay of the profile. In view of the estimate just mentioned and arguing as in 
Lemma[5l we have < C/(l + l^;]^)- Therefore, using the multipole expansion in 

(|72p and the decay property of G{x) from Lemma [A. H we conclude that 

p^9,^j^+ o(|x|-2), a = , / f{y) dy. (74) 



27r cos2 Oh Jr 

Clearly, a > in (l74l) . and to complete the proof we need to show that a > 0. To 
establish this fact, we first note that fi{x)dx > and Jjj f2{x)dx = 0. Indeed, it is 
easy to see that fi{x) > for all a; € R, and the second equality follows from the fact 
that the operator {—d'^/dx'^Y/'^ is self-adjoint and that constants lie in its kernel. To 
show that /jjj f^{x)dx > 0, we use ([T]) and symmetrize the integral to obtain 

1^ f f {sin p{x) -sin p{y)) {cos p{x) -cos p{y)) 
i^{x)dx ^ -— I I dxdy 



^ r ^sin(p(:r)+p(y))sin2(£M_PM 
27r {x - yf 



dx dy, (75) 



where we used trigonometric identities to arrive at the last line. In view of the fact 
that p e {9h: ^tt), the right-hand side of ([75l) is positive, and the claim follows. 

This concludes the proof of the theorem. □ 

Let us remark that the arguments in the proof of uniqueness above, with 
the test function d*'{x) := arcsin(i sin'!9'^^^(x) + (1 — i) sin^?(^^(a;)) for x > and 
1?* = TT - arcsin(tsini?(i)(a;) + (1 - t) sini?(^)(a;)) for a; < 0, where t G [0, 1], could 
also be used to prove uniqueness of the critical points of the energy E in the class 
of solutions of pTjl with values in (??h,7r — 'dh), obeying and crossing the value 
of ^ only once at the origin (for a similar treatment see [M]). However, since the 
computations in this case become exceedingly tedious and the precise behavior of the 
solutions at infinity may be needed, we have not pursued this question any further in 
this paper. Nevertheless, establishing such a uniqueness result would be helpful for 
interpreting the results of the numerical solution of ([TT|) as the Neel wall profiles. 

It would also be interesting to see if the one-dimensional Neel wall profiles 
considered in this paper are the only minimizers (or even critical points) of the two- 
dimensional thin film energy in Q with respect to perturbations with compact support 
that have the asymptotic behavior given by ([12]). We note that in the case v = this 
problem reduces to the famous problem of De Giorgi, whose solution in two space 
dimensions is now well understood [3S] (see also [351 for a recent overview) . Whether 
such a result remains valid in the presence of a non-local term {v > 0) remains to 
be seen (one result in this direction was obtained in [H]). Let us note that while 
for the local problem a continuous family of solutions obtained by rotations of the 
one-dimensional profile exists, in the non-local problem the orientation of the wall is 
expected to be fixed by the condition that the net charge across the wall be zero. The 
latter only allows walls that are parallel to the easy axis. 
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Appendix A. 

The following lemma establishes the basic properties of the fundamental solution for 
the operator L (see also [TH Sec. 5.1]). 

Lemma A.l. Let G{x) he the fundamental solution for the operator L defined in 
(pj) . Then 

G(x) = / ^ rrrdt. (A.l) 

In particular, G G C°°(R\{0}) n Lip(R) n L°°(R) n L^R); G > 0, G{x) = G{-x), and 

G(^) = ^H^l^r' + 0{\x\-'). (A.2) 
ZTT cos^ Oh 

Proof. The proof is a simple application of Fourier transform and contour integration 
techniques, which we present here for completeness. Observe first that the Fourier 
transform G{k) = Jjj e~'^^^G{x) dx of G{x) is well-defined and is given by 

^^^^ " |fc|2 + i,.COs2 eh\k\+C0s2 Oh ^^'^^ 

Interpreting |fc| = as an analytic function of k in the complex plane with a branch 
cut on the imaginary axis, i.e., defining •^/(a; + jy)^ = \x\ + iysgnx for x 7^ 0, wc can 
write the formula for inverting the Fourier transform of G as 

1 r ikx 

27r Jr + ijycos^ OhVk^ + cos^ Oh 

and treat it as a contour integral along the real axis. In particular, G{x) is even, and 
in the following it suffices to consider only a; > 0. 

It is easy to see that with our choice of the analytic branch the integrand in (|A.4p 
has no poles. Furthermore, since the integral in (IA.4I) over a semicircle of radius i? > 
in the positive imaginary half-plane vanishes for a; > as i? — >■ 00 , we can deform the 
contour of integration to run back and forth along the positive imaginary axis. After 
some algebra, we then find that the integral in (jA.4[) coincides with that in (jA.l[) . 

From the representation in (jA.ip . one immediately concludes that G{x) is positive, 
bounded, smooth for all a; 7^ and Lipschitz-continuous at a; = 0. The estimate 
in (jA.2p is then obtained by an elementary asymptotic analysis of the integral in 

dsn). □ 
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